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On the dynamics of the Meissner effect 
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The question of how a metal becoming superconducting expels a magnetic field is addressed. It 
is argued that the conventional theory of superconductivity has not answered this question despite 
its obvious importance. We argue that the growth of the superconducting into the normal region 
and associated expulsion of magnetic field from the superconducting region can only be understood 
if it is accompanied by motion of charge from the superconducting into the normal region. From a 
microscopic point of view it is shown that the perfect diamagnetism of superconductors requires that 
superconducting electrons reside in orbits of spatial extent 2An, with An the London penetration 
depth. Associated with this physics, the spin-orbit interaction of the electron magnetic moment 
and the positively charged ionic background gives rise to a “Spin Meissner” effect, the generation 
of a macroscopic spin current near the surface of superconductors. We point out that both the 
Meissner and the Spin Meissner effect can be understood dynamically under the assumption that 
the superfluid condensate wavefunction T (r) does not screen itself, just like the T(r) for an electron 
in a hydrogen atom. We argue that the conventional theory of superconductivity cannot explain 
the Meissner effect because it does not contain the physical elements discussed here. 

PACS numbers: 


I. INTRODUCTION 

It is generally [1] believed that the conventional theory 
of superconductivity explains the Meissner effect [2] , the 
telltale property of superconductors. This erroneous be¬ 
lief arises from two related widely held misconceptions. 
The first misconception is that the Meissner effect is an 
“equilibrium phenomenon” [3], that results from a par¬ 
ticular form for the relation between current density and 
magnetic vector potential in the equilibrium state of a 
superconductor [4-6]. The second misconception is that 
for a theory of superconductivity to predict the Meiss¬ 
ner effect it is sufficient that it yields a lower free energy 
for the system when the magnetic field is excluded than 
when the magnetic field remains inside the material [7]. 

The second misconception is easily dispelled by a coun¬ 
terexample. A superconducting ring with a finite current 
that generates an enclosed quantized magnetic flux has a 
higher free energy than the ring with no current and no 
magnetic flux. Yet there is no mechanism for the system 
to reach this lower energy state, and it will remain in the 
finite current “metastable” state forever. 

The first misconception is also easily disproved. The 
Meissner effect is the process by which a metal cooled 
below its critical temperature in the presence of a mag¬ 
netic field in its interior expels the magnetic field and 
reaches the equilibrium superconductive state with the 
magnetic field excluded. A calculation that deals only 
with quantities in the final state of this process, such as 
performed in Refs. [3-6] and many others cannot (by def¬ 
inition) prove that the system will reach this final state 
nor explain how the system does this. In the absence of a 
description of the process by which the equilibrium state 
is reached one can equally well conclude instead that the 
theory predicts that the system will forever remain in 
the initial “metastable” state, just like the ring discussed 
above. 


Thus we argue that the current theory of supercon¬ 
ductivity does not explain nor predict the Meissner ef¬ 
fect. Rather, the most natural conclusion one should 
infer from the theory is that a magnetic field either in 
the interior of a ring, a hollow sphere or a solid body 
will remain unchanged (except perhaps for tiny adjust¬ 
ments to account for quantized flux) when the system is 
cooled below its critical temperature. That this should 
happen for all three cases was “proven” theoretically by 
Lippmann [8] in 1919 based on Faraday’s law. It was 
confirmed experimentally by Onnes and Tuyn [9] in 1924 
for the case of a hollow sphere. But it was disproved 
experimentally by Meissner and Ochsenfeld in 1933 [10] 
for the case of a solid body. We argue that how super¬ 
conductors are able to prove Lippmann wrong and ignore 
Faraday’s law is not explained by the conventional theory 
of superconductivity. 

There has been a recent claim in the literature that the 
Meissner effect needs no explanation because any per¬ 
fect conductor will show a Meissner effect [11]. echoeing 
a similar claim made earlier [12]. These claims are in¬ 
correct, as has been proven by several authors recently 
[13, 14] as well as earlier [15, 16]. For example, Ref. [14] 
shows that for the Meissner state to be energetically fa¬ 
vorable requires the lowering of energy achieved by the 
superconducting condensation energy, which would not 
exist for a perfect conductor. While this certainly proves 
Ref. [11] ’s claim wrong it does not prove that the Meiss¬ 
ner effect is explained by the conventional theory for the 
reasons given above. 

Why is the question whether or not the conventional 
theory of superconductivity explains the Meissner effect 
relevant? Because if the conventional theory cannot ex¬ 
plain the Meissner effect it cannot be the correct the¬ 
ory of superconductivity, despite its many apparent suc¬ 
cesses. We argue that the conventional theory lacks es¬ 
sential physical ingredients that are necessary to explain 
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the Meissner effect. 

After discussing in more detail the absence of an expla¬ 
nation of the Meissner effect in the conventional theory 
in the next section, in the remainder of this paper we 
discuss what is required of a theory of superconductivity 
that can explain the Meissner effect. The theory of hole 
superconductivity contains those physical elements [17]. 


II. CONVENTIONAL UNDERSTANDING OF 
THE MEISSNER EFFECT 


The London equation, that provides a phenomenolog¬ 
ical description of the Meissner effect, is usually made 
plausible [18-20] by starting from the ‘acceleration equa¬ 
tion’ for a perfectly conducting fluid 


dv s _ q - 
dt m 


( 1 ) 


with v s the superfluid velocity for carriers of mass m and 
charge q and E the electric field, using Faraday’s law 


d_ 

dt 


(V x v s ) 


q dB 
me dt ’ 


( 2 ) 


integrating Eq. (2) in time and setting the integration 
constant equal to zero to yield the London equation 


V x v a = ——B 
me 


( 3 ) 


with B the magnetic field. Combined with Maxwell’s 
equation V x B = (4,it/c)J , with J = n s qv s the current 
density, Eq. (3) yields 


V 2 B = 


4irn s q 2 

me 2 




( 4 ) 


describing the fact that a magnetic field cannot penetrate 
a superconductor beyond a distance A l from the surface. 

However, integrating Eq. (2) in time yields instead of 
Eq. (3) 


VxiT s (r, t)—Vxi7 s (r, t = 0) = — —<B(f,t)—B(f,t = 0)). 

me 

( 5 ) 

If a normal metal is cooled into the superconducting 
state in the presence of a spatially uniform magnetic field 
B(r,t = 0) = Bq throughout its interior, the initial su¬ 
perfluid velocity v s (f,t = 0) = 0 and Eq. (5) yields 


V x v s (r,t) = —— (B(f, t) — Bo) (6) 

me 

which is not the same as Eq. (3). Quite the contrary, 
Eq. (6) implies that v s {r,t) = 0 and B(r,t) = B 0 for 
all times i > 0, so that the magnetic field remains un¬ 
changed inside the superconductor. Thus, this derivation 
certainly does not describe a process by which the sys¬ 
tem will achieve the Meissner state described by Eq. (3), 
rather it predicts that such a state will never be achieved. 


Within the conventional (BCS) theory of superconduc¬ 
tivity the Meissner effect is ‘proven’ [4-6] by calculating 
the London Kernel K ( q) relating the Fourier components 
of the current density J(q) and a static magnetic vector 
potential A(q) 


J(q) = ~^ k (q)Aq) O) 

and showing that 

K(q -A 0) = ± ^ 0 (7b) 

when the system is described by the BCS wavefunction. 
Eq. (7) is equivalent to Eq. (3). A great deal of literature 
was generated around the question whether Eq. (7) could 
be proven in a gauge-invariant fashion, and eventually 
this was achieved to everyone’s satisfaction [21]. However 
as discussed in the introduction this does not address the 
key question which is, how does the system achieve the 
BCS state that satisfies Eq. (7) starting from the initial 
state that doesn’t? 

Furthermore, there is an inherent problem in the lin¬ 
ear response argument Eq. (7). It implicitly assumes 
that the system is initially in the BCS state, then a uni¬ 
form magnetic field is created in its interior, and sub¬ 
sequently the system responds to it by generating the 
Meissner current that cancels the interior magnetic field. 
However it is physically impossible to ‘create’ a uniform 
magnetic field inside a material without magnetic lines 
‘cutting through’ the material, since Maxwell’s equations 
do not allow magnetic field lines to be created out of 
nowhere, they are always necessarily closed. If the sys¬ 
tem is initially in the BCS state it is a perfect conductor 
and magnetic field lines cannot cut through it so the lin¬ 
ear response situation cannot be set up. In addition, a 
system that is already in the BCS state cannot have a 
magnetic field in its interior because this is incompatible 
with global phase coherence. What actually happens in 
the Meissner effect is that the system in the normal state 
is cooled below T c , it is not in the BCS state initially, 
and in the process of entering the BCS state and estab¬ 
lishing macroscopic phase coherence the magnetic field is 
expelled. This complex process is not described by Eq. 

( 7 ). 

A rationale for the conventional argument seems to 
be [22] that when the system is cooled below T c in the 
presence of a weak uniform magnetic field, it is expected 
that thermal fluctuations will lead to the transition into 
the superconducting state with the field expelled because 
it has lower free energy than the normal state. However, 
we argue that thermal fluctuations are local and such 
random fluctuations cannot generate the global surface 
current necessary to expel the magnetic field except with 
vanishing probability for a macroscopic system. 

A somewhat more general ‘proof’ of the Meissner ef¬ 
fect [23] starts from the assumption that the many- 
electron superfluid condensate of a superconductor can 
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be described by a complex macroscopic wavefunction 
T(r). This was first done in a phenomenological way by 
Ginzburg and Landau [24], later it was shown that such 
a description can be derived from BCS theory under cer¬ 
tain conditions [25, 26], and the predictions of Joseph- 
son [27] and subsequent experimental verifications [28] 
established without doubt that it is a correct description 
of the superfluid condensate that captures the essence of 
superconductivity, whether BCS theory is valid or not. 
The amplitude of this macroscopic wave function T (r) is 
related to the density of superconducting carriers n s 

$(7) = |*(r}|e^ = nJ/V**) (8) 

and the gradient of the phase cj>(r) is related to the su¬ 
perfluid velocity v s according to the relation 

v s = —V<^>- —A. (9) 

m me 


within BCS-Bogoliubov-de-Gennes theory the supercon¬ 
ducting order parameter will not relax spontaneously to 
its equilibrium value. 

III. HOW A PERFECT CONDUCTOR EXPELS 
A MAGNETIC FIELD 

Let us consider again the equation of motion for a per¬ 
fectly conducting fluid. Eq. (1) is not quite correct be¬ 
cause it ignores the difference between total and partial 
time derivatives, and because the magnetic Lorentz force 
is omitted from the right side. Taking both of these facts 
into account leads instead of to Eq. (2) to the equation 
[35] 

dui 

— = v x x w) (10a) 

with 


The wavefunction 4'(r) describes the state of all the 
Cooper pairs in the system, that share a common phase 
</>(r) that is coherent over macroscopic distances [28]. 
Many of the unique and universal properties of super¬ 
conductors such as flux quantization and the variety of 
phenomena exhibited by Josephson junctions and ‘weak 
links’ follow from this simple macroscopic description. 
Upon taking the curl on both sides of Eq. (9) the Lon¬ 
don Eq. (3) results, hence this is assumed to be a proof 
of the Meissner effect. However once again, in assuming 
that the system is described by the macroscopic wave- 
function Eq. (8) one is assuming that the magnetic field 
has already been expelled, since Eq. (8) is not valid in the 
presence of a magnetic field in the interior of the super¬ 
conductor, and the process by which the system reaches 
or doesn’t reach this state starting from the initial state 
with the magnetic field in the interior is not discussed. 

Finally, an ‘energetic’ argument for the Meissner effect 
is given for example in ref. [7], where it is shown that 
the magnetic field distribution where the magnetic field 
is excluded except in a region within Al of the surface 
minimizes the free energy. However as already empha¬ 
sized, this does not explain the process of field expulsion 
nor predicts that the state will be reached. 

In summary, we argue that these arguments, which are 
generally assumed to prove that the conventional theory 
of superconductivity describes the Meissner effect, in fact 
do not do so, and leave completely open the question 
whether or not the conventional theory of superconduc¬ 
tivity can describe the Meissner effect. 

More recently, there have been calculations describing 
the normal-superconductor transition in a magnetic field 
[29, 30] using time-dependent Ginzburg Landau theory 
[31, 32]. In TDGL theory it is assumed that the su¬ 
perconducting order parameter (4/(r)) relaxes exponen¬ 
tially to its equilibrium value in a non-equilibrium sit¬ 
uation. However, no justification for this assumption is 
presented. In fact, very recent work claims [33, 34] that 
this assumption of TDGL theory is incorrect and that 


w = V x v s + —B (10b) 

me 

the ‘generalized vorticity’. The London Eq. (3) is 

w(r,t) = 0 (11) 


and the initial condition at the moment a system is cooled 
into the superconducting state in a uniform magnetic 
field Bq is 

w(f, t = 0) = w 0 = —Bo- (12) 

me 


To explain the Meissner effect along these lines one has 
to explain how w evolves in time from its initial value 
Eq. (12) to its final value Eq. (11) at all r’s at a later 
time. 

Let us assume for simplicity a cylindrical geometry 
with B in the £ direction, w = wz and Eq. (10a) can be 
rewritten as 


dw 

~dt 


Id 

— — [rwv r ) 
r or 


(13) 


which indicates that a radial velocity v r ^ 0 of the fluid 
is a prerequisite for w to change in time. In the absence 
of a radial velocity, dw/dt = 0 and w(f,t) = u>o for all 
times and the magnetic field in the interior of the metal 
becoming superconducting would remain unchanged, in 
agreement with Lippmann [8] and in contradiction with 
experiment. 

Eq. (13) can be rewritten as 


dw 

~dt 


-V • (wv s ) 


(14) 


which is a continuity equation. It says that for w to 
change locally it has to be carried away by the fluid flow. 
Integrating Eq. (14) over a cylinder of radius R and using 
Gauss’ theorem yields 



drr—w(r, t) 


—Rw(R 1 t)v r {R 1 1 ) 


(15) 
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and integrating over time and assuming at time T the 
Meissner state w(r, T) = 0 for all r < R has been reached 
yields the condition 

2 f T 

wo = — dtw(R,t)v r (R,t). (16) 

R J o 

According to this calculation in order to achieve the 
Meissner state in the interior of a cylindrical supercon¬ 
ductor of radius R , i.e. w(r) = 0 for all r < R, the entire 
superfluid has to flow out of the material carrying the 
generalized vorticity w with it. 

Of course superconductors don’t do that. The reason 
this calculation does not apply is that it assumes the sys¬ 
tem first becomes a perfect conductor and subsequently 
expels the magnetic field As we discuss in the next sec¬ 
tion, this is not what occurs in reality. Nevertheless this 
analysis is useful because it indicates that an outflow of 
charge is necessary to understand the dynamics of the 
Meissner effect. 


IV. HOW THE TRANSITION OCCURS 

In a recent paper [36] we have discussed the kinetics of 
the normal-superconductor transition in a magnetic field. 
Figure 1 shows three generic scenarios. As discussed in 
[36], the simplest scenario (a), where a current develops 
near the surface of the sample and the magnetic field is 
uniformly expelled can be excluded on theoretical as well 
as experimental grounds. In scenario (b) superconduct¬ 
ing regions nucleate at random and expand. Scenarios 
(b) and (c) are in essence similar so we will limit our¬ 
selves to scenario (c), where the superconducting phase 
expands from a nucleus at the center towards the surface 
with a surface current that excludes the magnetic field in 
the interior of the superconducting region. This has been 
implemented experimentally in a cylinder by applying a 
slightly smaller magnetic field in the central region [37]. 

Both in scenarios (b) and (c) a Faraday electric field E 
exists at and near the boundary of the domains because 
the magnetic field is changing as the domains expand. 
The direction of this electric field is opposite to the di¬ 
rection of the boundary current as shown in Fig. 1(c), 
to induce a current opposing the change in magnetic flux 
(Lenz’s law). Such a counterclockwise electric field also 
exists throughout the interior of the sample in (a) during 
the transition. 

Focusing on scenario (c), an explanation of the Meiss¬ 
ner effect has to explain how this current at the surface of 
the expanding superconducting region is generated in the 
first place, how it is maintained on the face of the Faraday 
electric field that tries to suppress it, how its momentum 
and kinetic energy is transferred radially outward as the 
superconducting phase expands, and how this growing 
momentum (and angular momentum) is compensated so 
as to not violate fundamental conservation laws. None 
of these questions is addressed by theoretical treatments 



FIG. 1: Three conceivable routes for the magnetic field ex¬ 
pulsion in a cylindrical superconductor. The dots represent 
magnetic field lines coming out of the paper. The arrows give 
the direction of the currents (I). In (a), a current starts devel¬ 
oping within A_l of the surface of the cylinder that gradually 
increases in magnitude, causing a gradual uniform decrease 
of the magnetic field in the interior. In (b), domains of loop 
currents start at various random locations that nullify the 
magnetic fields in their interior, that gradually expand their 
radii and coalesce and merge with neighboring domains. In 
(c), a single circular domain starts at the center and expands 
in radius until it reaches the boundary of the cylinder. 


of this process based on the conventional BCS-London- 
TDGL theory of superconductivity [29, 30]. 

Figure 2 (a) shows more stages of the process of mag¬ 
netic field expulsion (c) of Fig. 1. BCS-London-TDGL 
theory would say that the physics driving the expansion 
of the field-free region is the superconducting conden¬ 
sation energy. The energy of the system is lowered as 
Cooper pairs form and become phase coherent, condens¬ 
ing into the macroscopic superconducting state described 
by ^(r). Because the establishment of phase coherence 
requires that no magnetic field exists in that region, the 
region has to become field free, and this occurs through 
a BCS order parameter spontaneously growing and re¬ 
laxing towards its equilibrium value [29, 30]. This is a 
phenomenological treatment that does not explain how 
the condensation into the BCS state is connected to the 
azimuthal motion of the carriers that needs to be gener¬ 
ated at the boundary, and in particular what is the force 
that drives these carriers to move in direction opposite 
to the electric force exerted on them by the Faraday elec¬ 
tric field. Within BCS-London-TDGL theory there is no 
radial charge motion associated with the outward mo¬ 
tion of the phase boundary between superconducting and 
normal regions. 

The behavior shown in Fig. 2 (a) would also occur if 
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(a) Growth of superconducting phase 



(b) Outward flow of perfectly conducting fluid 



FIG. 2: (a) More stages of the process of magnetic field 

expulsion as a system goes superconducting according to the 
scenario (c) of Fig. 1. At and near the boundary of the 
expanding field-free region there is an electric field pointing 
counterclockwise, due to the changing magnetic flux and Fara¬ 
day’s law, and a current flowing clockwise that cancels the ap¬ 
plied magnetic field in the interior, (b) Schematic depiction of 
a perfectly conducting fluid (o = oo) that flows radially out¬ 
ward with radial velocity v r . The carriers at the boundary 
experience a Lorentz force. Assuming the sign of the charge 
q is positive for definiteness the Lorentz force Fl = ( q/c)v r B 
points in the clockwise direction. The resulting electric cur¬ 
rent I at the boundary flows clockwise (for negative charge 
carriers the Lorentz force would be in opposite direction, the 
current in the same direction), generating a magnetic field 
opposite to the external field so that no magnetic field pene¬ 
trates the perfect conductor. 



I ' ’El Bo' 

iSfii 


H 


Xp(t) 




(b) 


Fe_ Fb ] F b 

noving superfluid 


FIG. 3: Magnetic field Bo points out of the paper, (a) 
Superconducting-normal phase boundary xo{t) moving in the 
+* direction with no associated charge flow in the x direction, 
(b) Superfluid moving in the +x direction with boundary at 
xo (t). For both cases the magnetic field is expelled from the 
region x < xo(t) by the current J y flowing within A l of the 
boundary, and a Faraday field E y exists due to the changing 
magnetic field. The electric Lorentz force (Fe) drives the car¬ 
riers in direction opposite to their direction of motion. Only 
for case (b) the carriers also experience a magnetic Lorentz 
force Fb driving them in the direction of their motion, oppo¬ 
site to Fe and of equal magnitude. 


a core of high density perfectly conducting material (in¬ 
finite conductivity) would expand radially outward, as 
shown schematically in Fig. 2 (b). Here the charge carri¬ 
ers are moving radially out and there is a Lorentz force Fl 
acting on them giving them an azimuthal velocity (in the 
clockwise direction assuming positive carriers) such that 
they generate a magnetic field opposite to the applied 
one, thus not allowing the magnetic field to penetrate 
the interior of the perfectly conducting region (except to 
within a distance Al of its surface). We discuss the pro¬ 
cess quantitatively in what follows in a planar geometry 
for simplicity. 


e v (x) = ^r^ e(x-xo)/Ai ( 17b ) 

with Bqz the applied magnetic field and E y in the pos¬ 
itive y direction. Assuming current carriers of charge q , 
mass m and density n s , their speed in the y direction is 

v y (x) = - C B 0 e ix ~ Xo) ^ L (17c) 

4irn s q\L 

parallel (antiparallel) to the current, which flows in the 
—y direction, for q > 0 (q < 0). For x > Xq we assume 
for simplicity 

B(x) = B 0 (18a) 


V. MEISSNER EFFECT IN PLANAR 
GEOMETRY 


Ey( X) 


B 0 dx 0 
c dt 


(18b) 


Figure 3 shows the processes of Fig. 2 in a planar 
geometry. The magnetic field points in the z direction, 
the phase boundary advances in the +x direction, and the 
Faraday electric field and Meissner current point in the y 
direction, xq (t) is the boundary of the superconducting 
region, moving upward at rate dxq/dt. The magnetic and 
electric fields for x < xq are given by 


(x—Xq)/\l ( 17a ) 


This assumption corresponds to the treatment of ref. [36] 
for the particular case where the magnetic field that is 
being expelled is very close to the critical field (p —> 0 in 
the notation of ref. [36]). 

Equations (17a,b,c) follow from the relation between 
current and velocity J = n s qv , the expression for the 
London penetration depth Eq. (4), Maxwell’s equations 


V x B = 


47T 


dB 


47T 


-j -> ,, Jy 

OX c 


4:Trn s q 


(19a) 


B(x) = B 0 e 


c 
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- - ldB dE v 1 dB 

V xE= —— ==> = — — 

c at ox c ot 

and London’s equation 

v x j=-™1b==>^ = ~J-b. 

me ox me 

The electromagnetic force acting on the carriers is 

dv q -> q _ -< 

to— = — E H- v x B 

at m me 

and its component in the y direction yields 


dv v 


dv v 


= -E„ - 


-v x B 


(19b) 


(19c) 


( 20 ) 


( 21 ) 


dt dx m " me 

For the case of Fig. 3(b) (moving superfluid) we have 

dx o 


v x = 


dt 


( 22 ) 


and the terms on the right-hand side of Eq. (21) satisfy 


Fe = -Ey = —v x B = F b (23) 

m toc 

so that the electric and magnetic forces in Eq. (21) ex¬ 
actly cancel out. The left-hand side of Eq. (21) is also 
identically zero from Eqs. (17c) and (22). Carriers are 
accelerated by the magnetic Lorentz force and deceler¬ 
ated by the electric force arising from Faraday’s field, 
which exactly cancel in steady state. Thus, a dynamical 
explanation of flux expulsion is provided by these equa¬ 
tions: the magnetic Lorentz force acting on the outflowing 
charge drives the current that nullifies the magnetic field 
in the interior against the Faraday electric force. 

For the processes in Figs. 3 to happen requires a driv¬ 
ing force in the +x direction. The way this works is very 
clear for case (b), moving superfluid: there is a magnetic 
Lorentz force on the J y current carriers in the —x direc¬ 
tion: 

F x = ?-{vx B) X = \B (24) 

c c 

so that an equal and opposite force needs to be applied in 
the +x direction for the superfluid to advance in the +x 
direction. This requires expenditure of energy per unit 
time per carrier —F x v x , so that the energy per unit area 
per unit time spent is, from Eqs. (17a), (17c) and (24) 


dxn s ( Fj^Vyfj 



B(x) 2 Bq dx o 

dx --— v x = -— 

47tAl 87t dt 


(25) 


which equals the rate of change of magnetic energy per 
unit area as the boundary moves, as required by energy 
conservation. There is no energy dissipated in this pro¬ 
cess. 

The same rate of energy expenditure is required for the 
process of Figure 3(a) to expel the magnetic field. This 


energy obviously is supplied by the superconducting con¬ 
densation energy. However since in this case v x = 0, Eq. 
(25) cannot be used to understand how the condensation 
energy causes the phase boundary to advance. The force 
F x given by Eq. (24) still exists in this case, but since 
carriers don’t move in the x direction this force does not 
do any work. The work done by the condensation process 
in moving the phase boundary in the -\-x direction with¬ 
out displacing charge carriers in the -\-x direction occurs 
through some unknown way. 

In addition since charge does not move in the direc¬ 
tion of the moving phase boundary there is no magnetic 
Lorentz force driving the Meissner current in Fig. 3(a) 
and Eq. (21) is not satisfied. Eq. (21) becomes 


dv v 9 F 
dt to v 


(26) 


which certainly does not describe Fig. 3(a) since it says 
that the current should be flowing in the opposite direc¬ 
tion, parallel rather than antiparallel to E y . In steady 
state the left side in Eq. (26) is zero and there is an 
uncompensated right-hand-side. 

One would have to include other dynamical effects to 
‘fix’ Eq. (21) and to give a dynamical explanation of the 
energy transfer process for the case of Fig. 3(a) other 
than Eq. (25). This has not been done within the con¬ 
ventional theory. We argue that the dynamics described 
by Eqs. (20) and (25), requiring charge motion in the 
same direction as the phase boundary motion, is the only 
physical way to describe the flux expulsion. 


VI. BACKFLOW 

As discussed in the previous section, the strong simi¬ 
larity in the physics shown in Figs. 3(a) and 3(b) is com¬ 
pelling evidence that in the process of a metal becoming 
superconducting an outflow of charge in direction nor¬ 
mal to the superconductor-normal metal phase boundary 
into the normal region takes place. This explains both 
how the carriers in the supercurrent flowing parallel to 
the phase boundary can move in direction opposite to 
the Faraday electric force, and how the superconducting 
condensation energy is used to do the work required to 
expel the magnetic field: the condensation results in an 
electromotive force [38] driving the carriers to move in 
the +x direction, i.e. from the superconducting into the 
normal region, opposing the electromagnetic force Eq. 
24. 

In order for this not to create an enormous charge im¬ 
balance it requires a backflow of normal carriers in the — x 
direction. Figure 4 shows schematically how this can hap¬ 
pen. The layer of superfluid of thickness A l next to the 
phase boundary with the velocity pattern given by Eq. 
(17c) moves forward a distance and an equal amount 
of normal fluid moves backward, and in the process be¬ 
comes superfluid. This satisfies energy and momentum 
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FIG. 4: Schematic depiction of how the superfluid-normal 
phase boundary moves into the normal region. Superfluid 
layer of thickness A l carrying the screening current J y moves 
forward and an equal amount of normal fluid moves backward 
and becomes superconducting. 



FIG. 5: As carriers become superconducting (s carriers) they 
thrust forward into the normal region over a boundary layer of 
thickness A l , and are deflected by the Lorentz force acquiring 
speed v y = — c/(Airn s q\L,)Ba. This process creates an electric 
field E x in the —x direction that drives normal carrier (n 
carrier) backflow. Here, “n” stands both for “normal” and 
“negative”. The normal carriers do not acquire a large v y 
in opposite direction because they scatter off impurities and 
transfer their y —momentum to the lattice. 


conservation and gives a simple description of the pro¬ 
cess. The driving force for the backflow is the electric 
field in direction perpendicular to the phase boundary 
created by the forward flow. 

However, considered in detail the process is likely to 
be more complicated. We propose that the forward mo¬ 
tion occurs for carriers that are in the process of be¬ 
coming superconducting, which are initially (essentially) 
not moving in the y direction and acquire a large veloc¬ 
ity parallel to the phase boundary through the magnetic 
Lorentz force. As discussed in Ref. [36], in order for a 
carrier moving in the +x direction to acquire the trans¬ 
verse velocity Eq. (17c) through the action of the mag¬ 
netic Lorentz force requires a displacement over a dis¬ 
tance A l in the x direction. Thus, we envision a process 
of forward flow and backflow over a boundary layer of 
thickness A^ as the phase boundary advances, as shown 
schematically in Fig. 5. As carriers become supercon¬ 
ducting they thrust forward a distance A l , which creates 
an electric field in the +x direction (assuming these are 
negatively charged carriers) causing backflow of normal 
carriers in this layer. As the phase boundary advances 
at rate vq = dxo/dt , normal (negatively charged) carri¬ 
ers in the boundary layer are back-flowing at speed Vo, 
or equivalently positive normal carriers (holes) move for¬ 
ward together with the phase boundary. Because nor¬ 
mal carriers scatter off the lattice they do not acquire a 




FIG. 6: Expansion of superconducting phase in a cylindrical 
geometry. On the left, the dark grey circle is superconducting 
and the light grey ring (of thickness ~ Az,) is about to become 
superconducting. As the ring becomes superconducting the 
current direction in the ring switches from the +9 to the —9 
direction (red arrows) which is opposite to the electric field 
direction. 


large v y from the action of the magnetic Lorentz force, 
instead they transfer their y —momentum to the lattice 
as a whole, thus accounting for momentum conservation. 

This gives a phenomenological description of the mag¬ 
netic field expulsion process that resolves the puzzles of 
the Meissner effect regarding forces and momentum con¬ 
servation. A complete microscopic description has not 
yet been proposed. A semiclassical explanation based 
on the proposal that superconducting carriers reside in 
mesoscopic orbits of radius 2A l [39] is discussed in the 
following sections. 

For completeness we show in figure 6 the current and 
fields in a cylindrical geometry. The magnetic field that 
is being expelled is H c (l—p) with p > 0, and it increases 
to H c at the phase boundary due to the magnetic field 
generated by the eddy current in the normal region. The 
equations are given in Ref. [36]. Again, the process is 
simply understood assuming the outward expansion of 
the phase boundary is accompanied by radial outflow of 
charge, driven by a radial electromotive force that pushes 
the charges against a radially inward pointing magnetic 
Lorentz force. An inward backflow of normal current 
(not shown) compensates for the charge imbalance. In 
the absence of radial outflow of charge there is no mech¬ 
anism to explain how the current reverses its azimuthal 
direction (red arrows) and flows against the direction of 
the electric field when the phase boundary advances, nor 
to explain how the growing angular momentum of the 
Meissner current is compensated [36, 42], 


VII. MEISSNER EFFECT AND 2A L ORBITS 

Consider a superconducting long cylinder in an applied 
magnetic field H. The Meissner current that nullifies the 
magnetic field in the interior flows within a distance A l 





















of the surface with speed 


v s 


— H. 

m e c 


(27) 


with e and m e the electron charge and mass. This is 
easily seen from the relation between magnetic field and 
vector potential for the cylinder, A = HAl- The Meiss¬ 
ner current density is 

71 6 2 

j s = n s ev s = ——A l H (28) 

m e c 


and from Maxwell’s equation V x H = (4n/c)j s it follows 
that 


n 4w ■ \ 
H = — j s Al 


c 

and combining Eqs. (28) and (29) yields 

1 47rn s e 2 
A| m e c 2 

The magnetization per unit volume is given by 




en s 


A lv s 


(29) 


(30) 


( 31 ) 


It is reasonable to assume that the magnetization results 
from a superposition of elementary magnetic moments jl 
resulting from the orbital motion of each superconducting 
electron: 


M = n s fl. 


(32) 


This yields for the magnetic moments 

eA l 

ft = - v s 

c 


(33) 


and from the general principle that the magnetic moment 
and orbital angular momentum £ of each electron are 
related by 


P = 


e t* 


2 m e c 


(34) 


we deduce that the angular momentum associated with 
each of these electrons is 

t- = 2meC ^ = m e v s ( 2X l ). (35) 

e 


as follows. For an electron in an orbit of radius r, ap¬ 
plication of an external field H yields through Faraday’s 
law a tangential electric field 


E = — 


r dH 
2c dt 


(36) 


and from the equation of motion m e dv/dt = eE, 


dv er dH 
dt 2 m e c dt 


(37) 


so that Eq. (27) results upon integration for r = 2A l- 
However there is a subtle question of self-consistency that 
we discuss in the next section. 


VIII. MEISSNER EFFECT AND MAGNETIC 
SUSCEPTIBILITY 

An external magnetic field H applied to a material 
gives rise to a magnetization M = \H and a total mag¬ 
netic field 


B = H + 4ttM = (1 + 4t t X )H. (38) 


A perfect diamagnet does not allow any magnetic field in 
its interior, hence is defined by y = — 1/(47t) according to 
Eq. (38). This is usually assumed to be the susceptibility 
of superconductors. 

However, Pippard points out [40] that this is a miscon¬ 
ception when applied to superconductors because “the 
mean field responsible for magnetizing an extended unit 
is not H but H”, from which it follows that M = yU 
rather than M = in Eq. (38), hence 


B = 


1 

1 - 47TX 


(39) 


and “perfect diamagnetism demands that y be infinitely 
negative” [40]. The same argument is made by Tinkham 

[ 41 ]- 

Consider a system of electrons of density n s per unit 
volume in orbits perpendicular to an applied magnetic 
field H. The Larmor diamagnetic susceptibility is given 

by 


n s e 2 
4m e c 2 


< 


r 2 > 


(40) 


where \/< r 2 > = r gives the spatial extent of the orbit, 
i.e. its ‘radius’. For 


Eq. (35) implies that the Meissner current that nullifies 
the magnetic field in the interior of a superconductor re¬ 
sults from each superfluid carrier moving with speed v s 
in a mesoscopic orbit of radius 2 Al. Just like in a mag¬ 
netic material, the superposition of elementary currents 
results in a macroscopic surface current. 

It can be seen that the fact that electrons reside in 2A^ 
orbits by itself is sufficient to give perfect diamagnetism, 


f=2A L (41) 

the diamagnetic susceptibility Eq. (40) is y = — 1/An 
according to Eq. (30), and we have argued in the past 
that this should apply to superconductors [42], Instead, 
Pippard and Tinkham argue, as discussed above, that 
f = oo is required in superconductors to give y = — oo 
and B = 0 (Eq. (39)). 
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Indeed it would appear at first sight that the Pippard- 
Tinkham point of view is the correct one. The orbits of 
radius given by Eq. (41) are highly overlapping since 


/ - 1/3 

2Al = nj 1/3 y —->> n~ 1/3 (42) 

V nr c 

where r c = e 2 /m e c 2 is the classical electron radius, much 
smaller than the inter-electron distance n s ■ Therefore 
one would expect the magnetic field generated by one 
orbit to reduce the magnetic field affecting nearby over¬ 
lapping orbits, and the net magnetic field resulting from 
application of an external magnetic field H to be, from 
Eq. (39) 


B = 


1 

1 — 47t(—1/47t) 



(43) 


which is not the complete flux expulsion observed in su¬ 
perconductors. 

The same result is obtained from the analysis in the 
previous section. Faraday’s law Eq. (36) should involve 
the total field B rather than H , hence instead of Eq. (37) 
we have 


dv er dB 

dt 2m e c dt 


(44) 


The total magnetic field is weaker than the applied mag¬ 
netic field due to the counterfield generated by the other 
electrons in the system: 

27rn p.v 

B = H + 47 tM = H + 47t rig/i = H - -—v (45) 

c 

using that y = —erv/2c for electrons with speed v in 
orbits of radius r. Replacing in Eq. (44) and performing 
the time integration yields 


v = — 


er 

2 m e c 


i + 


7T n s e 2 
m e c 2 


r 2 


H 


which for r = 2 Xl yields 


v = — 


eX l_ h _ Vs_ 
2 m e c 2 


(46) 


(47) 


predicting that electrons aquire only half of the required 
speed Eq. (27) to nullify the applied magnetic field H, 
hence that the applied magnetic field will not be nulli¬ 
fied in the interior but only reduced by a factor of 2, in 
agreement with Eq. (43). 

This result is puzzling because a perfect classical con¬ 
ductor will perfectly screen an applied magnetic field in 
its interior. We conclude that to make this semiclassical 
model agree with the classical behavior requires that Eq. 
(37) rather than Eq. (44) applies, in other words that 
electrons in the superfluid are only affected by the exter¬ 
nal magnetic field and not by the magnetic field generated 
by other electrons in the superfluid. 


Thus, to understand the perfect diamagnetism of su¬ 
perconductors in terms of 2A l orbits as suggested by Eqs. 
(32) and (35) it is necessary to assume that the mag¬ 
netic field created by superfluid electrons does not affect 
the superfluid itself. This is not an implausible assump¬ 
tion. It will occur if the superconducting electrons are 
described by a macroscopic wavefunction (f) that does 
not screen itself. Thus, when an external magnetic field is 
applied all the components of the macroscopic wavefunc¬ 
tion are subject to the entire applied field rather than to 
the external field reduced by the action of other compo¬ 
nents of the wavefunction. Under those conditions, Eq. 
(37) rather than Eq. (44) applies and the induced mag¬ 
netic field for electrons in orbits of radius 2Az, is precisely 
of the magnitude needed to completely cancel the interior 
magnetic field. 

This key property of the quantum-mechanical wave 
function was recognized from the outset of the develop¬ 
ment of quantum mechanics. Schrodinger expressed it 
clearly in his 1952 paper on the meaning of wave me¬ 
chanics [43]: “The original wave mechanical model of 
the hydrogen atom is not self-consistent. The electronic 
cloud effectively shields the nuclear charge towards out¬ 
side, making up a neutral whole, but is inefficient inside; 
in computing its structure its own field that it will pro¬ 
duce must not be taken into account, only the field of 
the nucleus. ” The reason Schrodinger spelled this out in 
detail is because it contradicted his intuition and physi¬ 
cal expectation that |tR(^| 2 would represent the charge 
density of the electron at position r rather than the 
probability of finding the electron at r. 

As we discuss in the next section, the same physics 
explains how a magnetic field is expelled from the interior 
of a system becoming superconducting, and furthermore 
we argue that this process cannot be explained in the 
absence of this physics. 


IX. DYNAMICAL EXPLANATION OF THE 
MEISSNER EFFECT 

The Landau diamagnetic susceptibility of conduction 
electrons is given by 

X = -^Bg(e F ) (48) 

with g(e F ) the density of states at the Fermi energy and 
ys the Bohr magneton. Using the expression for the 
free electron density of states g{e F ) = 2>n/2e F it is found 
that the Landau susceptibility Eq. (48) is the Larmor 
diamagnetic susceptibility Eq. (40) for 

< r 2 >= kp 2 . (49) 

and n = n s , hence that in the normal state electrons re¬ 
side in microscopic non-overlapping orbits of radius kp 1 . 
Therefore, in the process of becoming superconducting 
and establishing phase coherence electrons expand their 
orbits from radius kp 1 to radius 2Xl- 
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FIG. 7: A single charge carrier of positive charge expanding 
its orbit in a magnetic field pointing out of the paper acquires 
an azimuthal velocity in the clockwise direction, and itself 
generates a (small) magnetic field pointing into the paper. 



FIG. 8: Left panel: In the superconducting region the carriers 
reside in orbits of radius 2A l , in the normal region in orbits of 
radius kp 1 . At the boundary of the normal-superconducting 
region (dotted circle) the orbits expand from radius kp 1 to 
radius 2A l, causing magnetic field expulsion. Right panel: 
in the superconducting region the charge distribution is not 
homogeneous, there is an excess negative charge within A; 
of the surface of the superconducting region that spills over 
into the normal region, and an excess positive charge in its 
interior. 


The process of expansion of electronic orbits from ra¬ 
dius kp 1 to radius 2 Al provides a dynamical explanation 
of the Meissner effect. Due to the action of the Lorentz 
force on the radially outmoving electron it acquires an 
azimuthal velocity given by [44] 


v g = - 


qr 
2 me 


H 


(50) 


which is identical to Eq. (37), the speed acquired by a 
charge in an orbit of radius r when the applied magnetic 
field is increased from 0 to its final value H. The reason 
is Faraday’s law for Eq. (37), Lorentz force for Eq. (50). 
When the radial motion is over a distance r = 2Al the 
azimuthal velocity acquired is 


vg =- H 

me 


(51) 


which is the same as the speed of the Cooper pairs in 
the Meissner current Eq. (27). Note that it is important 
in this analysis that the magnetic field imparting the az¬ 
imuthal velocity to the electron in the expanding orbit is 
H rather than B. 

Figure 7 shows schematically the expansion of a single 
orbit in a magnetic field pointing out of the paper. As 
the orbit expands the carrier’s orbit cuts through mag¬ 
netic field lines and in so doing acquires an azimuthal 
velocity, shown clockwise in Fig. 7 assuming the carrier 
has positive charge. When the radius reaches 2 Al, the 
azimuthal speed reaches the value Eq. (27). In turn, the 
motion of this charge generates a magnetic field in the 
direction opposite to the applied field, i.e. into the paper 
in Fig. 7. For a single carrier of course the magnitude of 
this counterfield is negligible. 

Now we consider the growth of the superconducting 
region in a cylinder. In the normal outer region the 
orbits are microscopic, in the interior superconducting 


region the orbits have radius 2A l, and at the bound¬ 
ary between superconducting and normal regions the or¬ 
bits are expanding. This is shown schematically in Fig¬ 
ure 8 (left panel). The carrier in each expanding orbit 
reaches the final azimuthal velocity Eq. (27) when the 
radius of the orbit reaches its final value 2A^. Superpo¬ 
sition of these motions results in a current being carried 
within a layer of thickness A l from the boundary of the 
normal-superconducting region. As the superconducting 
region expands, this boundary current expands with the 
boundary, and no net current remains in the interior re¬ 
gion due to cancellation of the internal motions. When 
the boundary reaches the boundary of the sample, the 
system reaches the superconducting state where all the 
magnetic field has been excluded except within a layer 
of depth Al from the surface, and a Meissner current 
flows in this layer. In this process it is crucial that each 
carrier is affected by the full external magnetic field as 
its orbit expands, rather than one that is partially com¬ 
pensated by the magnetic field created by the expanding 
orbits of other carriers, in order that its azimuthal speed 
reaches the final value Eq. (51). Note in particular that 
within the enlarged orbits at the phase boundary in Fig. 
8 (left panel) there are small normal orbits that will ex¬ 
pand next, in a net magnetic field that is already weaker 
because of the field generated by the larger orbits enclos¬ 
ing them. However these smaller expanding orbits should 
not be affected by the magnetic field (in direction oppo¬ 
site to the applied field) created by the already enlarged 
orbits, because these carriers are all becoming part of the 
same 'F(r), hence they do not affect each other. 

To clarify this point further, we show in Fig. 9 the 
process of expansion of a single orbit. As the orbit cuts 
through the magnetic field lines the electron acquires 
azimuthal velocity vg due to the Lorentz force, which 
reaches the value Eq. (51) when the orbit reaches radius 
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ductor rotating with angular velocity wo is 


B = 


2m e c _ 

-w 0 

e 


(52) 


The resulting magnetization for a cylinder rotating 
around its axis is 



m e c 

2ne 


w o 


(53) 


FIG. 9: A single orbit expanding in the magnetic field. As 
the orbit cuts through the field lines (left panel) it acquires 
azimuthal velocity ve reaching the value Eq. (27) when the 
radius reaches 2A l- When the magnetic field lines move out ( 
right panel) they do not affect the azimuthal velocity because 
different parts of \l/(f) do not affect each other. 


2Al. Then, the magnetic field lines move out cutting 
through the 2A l orbit because of the compensating mag¬ 
netic field generated by all the carriers in 'F(r). This 
would slow down the azimuthal motion due to Faraday’s 
law if this was a normal carrier. Instead, because the 
carrier is part of ’F (f) its azimuthal speed is not affected 
as shown in the right panel of Fig. 9. In reality, the 
processes shown on the left and right panels of Fig. 9 do 
not occur sequentially but simultaneously. 

We conclude from this analysis that in the transition to 
superconductivity, in addition to the Cooper pairing pre¬ 
dicted by the conventional theory, the individual carriers 
forming the Cooper pair expand their orbits from mi¬ 
croscopic dimension (kp 1 ) to radius 2A^, without being 
affected by the magnetic field generated by other carri¬ 
ers becoming part of the same macroscopic wavefunction, 
and that this physics resolves the puzzle of the Meissner 
effect. 

In the right panel of Fig. 8 we show schematically 
the charge distribution. Because the expansion of the 
orbits involves some outward charge motion, our theory 
predicts that associated with the 2 Al orbits there is an 
excess negative charge near the boundary of the super¬ 
conducting region that spills over into the normal region 
[45, 46]. Thus the outward motion of the phase bound¬ 
ary can be understood both as associated with outward 
motion of negative charge into the normal region and as¬ 
sociated with enlargement of orbits from radius kp 1 to 
radius 2A^. The driving force for expansion of the orbits 
and associated negative charge expulsion is lowering of 
quantum kinetic energy, which is what drives supercon¬ 
ductivity according to the theory of hole superconductiv¬ 
ity [47, 48], 


X. 2A L ORBITS AND THE LONDON MOMENT 

The London moment is the magnetic moment gener¬ 
ated by a rotating superconductor [35]. The magnetic 
field that exists throughout the interior of a supercon¬ 


Assuming as in Eq. (32) that each supercarrier con¬ 
tributes magnetic moment fl to the magnetization yields 


h = 


m e c 

2nen s 


— \ 2 
wo — —a l u; o 

c 


(54) 


and the relation Eq. (34) between magnetic moment and 
angular momentum yields 


£ = m e (2\ L uj 0 )(2\ L ) (55) 


which describes carriers in orbits of radius 2 Al and tan¬ 
gential velocity 2 AlWo- This provides additional support 
to the interpretation that carriers in the superconduct¬ 
ing state reside in orbits of radius 2 Al. The fact that the 
London moment is parallel rather than antiparallel to the 
angular velocity demonstrates that the carriers forming 
the superfluid have negative charge [49]. 


XI. SPIN MEISSNER EFFECT 

The Spin Meissner effect [39] is the spontaneous gen¬ 
eration of a spin current within a London penetration 
depth of the surface of a superconductor when a metal is 
cooled into the superconducting state, given by 

J a =n s v^ = -n s ^ a x n (56) 

t±m e AL 

where h is the outward-pointing normal to the surface 
of the superconductor and a is parallel to the surface. 
The magnitude of the mass transport associated with 
this current in each direction is half the mass transport 
of the critical current of the superconductor. When a 
magnetic field B is applied the spin current component 
with a parallel to B slows down and the one with opposite 
spin direction speeds up, and when the magnetic field is 
such that the slower spin current component is stopped 
the superconducting state is destroyed [39]. 

The negative charge expulsion predicted by our theory 
[45] has as consequence that an outward-pointing elec¬ 
tric field exists in the interior of superconductors at suffi¬ 
ciently low temperatures. Because of this, a spin current 
originating in the Rashba spin-orbit interaction is ex¬ 
pected. However, the magnitude of the spin current given 
by Eq. (56) is orders of magnitude larger than would be 
expected from the ordinary Rashba effect [46, 50]. 
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The spin-orbit interaction of an electron in the pres¬ 
ence of an electric field E obtained from the Dirac Hamil¬ 
tonian is 

= (57) 


This can be represented by the Aharonov-Casher vector 
potential [51, 52] A a in the single-particle Hamiltonian 


H = 

Xr 


(p- -X ) 2 

(58a) 

c 


h -» 

(58b) 

- a x E. 


4m e c 


The term linear in A„ from Eq. (58a) yields Eq. (57) 
(for an interpretation of the term quadratic in A a see 
ref. [53]). Just like the ordinary vector potential A gives 
rise to a magnetic field B = V x A, the spin-orbit vector 
potential A a gives rise to an effective magnetic field [39] 

B„ = V x Xt = T -— (y-E)a = — pa (59) 

4m e c m e c 

that imparts an azimuthal velocity to the carriers in 
the expanding orbits, just as the ordinary magnetic field 
does [44], Here, p = V • E/4ir is the charge density that 
gives rise to the electric field E with which the magnetic 
moment of the electron interacts. 

The question now is, what is this charge density pi The 
superfluid density n s has associated with it a charge den¬ 
sity qn s , with q the charge of an individual carrier. Ex¬ 
periments such as the London moment, the gyromagnetic 
effect and the Bernoulli effect indicate that the charge 
carriers in the superfluid are electrons, hence q = e [49]. 
Again we have to assume that the superconducting fluid 
described by 'l'(r) does not screen itself and as a con¬ 
sequence the spin-orbit interaction affecting each super¬ 
fluid carrier results from the electric field generated by 
the full compensating ionic charge density 


P = \e\n s 


(60) 


rather than from the slight net charge imbalance resulting 
from charge expulsion, which is smaller than Eq. (60) by 
a factor u°/c [46]. Replacing Eq. (60) in Eq. (59) and 
using Eq. (30) for Al yields 

Trfr , , he _ . . 

B ° = = -saT (61) 

When the orbit expands to radius 2Al, the azimuthal 
speed acquired is 


v o - cXl Bu - h 

m e c 4TO e Ai 


(62) 


and just like for the Meissner effect, the internal motions 
cancel out and a spin current remains within Al of the 


surface, given by Eq. (56). The direction of the spin 
current is as given in Eq. (56). The angular momentum 
of electrons in 2A l orbits with speed given by Eq. (62) 
is 

£ = m e v°x( 2A L ) = |, (63) 

the same as the intrisic electron angular momentum due 
to spin. 

The condition Eq. (63) presumably has a topologi¬ 
cal origin and is what determines that the orbits expand 
to radius 2A^, which coincidentally is precisely what is 
needed to generate a magnetic field of just the right mag¬ 
nitude to cancel the external magnetic field and give rise 
to the Meissner effect. The driving force for the orbit 
expansion is lowering of quantum kinetic energy in the 
transition to the superconducting state [47, 53]. We re¬ 
gard the fact that the result Eq. (63) results from this 
analysis to be compelling evidence in favor of the validity 
of this model for the description of real superconductors. 


XII. MACROSCOPIC PHASE COHERENCE 

Macroscopic phase coherence is a hallmark of super¬ 
conductivity [2, 23, 27, 28]. The BCS wavefunction ex¬ 
hibits macroscopic phase coherence, however it does not 
provide an intuitive picture of what macroscopic phase 
coherence means, nor how it is established in the transi¬ 
tion from the normal to the superconducting state, nor 
how it is robustly maintained over macroscopic distances 
in the superconducting state, nor how the establishment 
of phase coherence is related to the Meissner effect. 

Instead, the theory discussed here provides a uni¬ 
fied explanation for how superconductors expel magnetic 
fields and how macroscopic phase coherence is estab¬ 
lished. Within our theory the superfluid wavefunction 
is composed of paired orbits of spin up and spin down 
electrons, each orbit of radius 2 A l, with distance between 
the centers of the orbits £, the superconducting coherence 
length [39]. We can think of the “phase” as a point in 
the electron’s orbit that is rotating with angular velocity 
w = v®/(2Al) = \e\/(2m e c)B cr . As the carriers condense 
into the superconducting state their orbits expand and 
overlap with each other, and this gives rise to phase co¬ 
herence because overlapping orbits have to have the same 
phase to avoid collisions, as shown schematically in Fig. 
10. We can easily understand that this phase coherence 
has to extend over the entire region occupied by the su- 
pcrfluid wavefunction ’l'(r). 


XIII. SUMMARY AND DISCUSSION 

In this paper we have argued that the generally held 
view that the conventional theory of superconductivity 
describes, explains and predicts the Meissner effect is 
incorrect. The Meissner effect is the process by which 
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FIG. 10: The “phase” is depicted as a black circle on the 
orbit. In the normal (n) phase (small orbits of radius kp 1 ) 
there is no correlation between the phases of different orbits 
because they are non-overlapping. In the superconducting (s) 
phase (large orbits of radius 2A l) the orbits overlap strongly 
and the phases have to be the same in the different orbits 
to avoid collisions. As time (t) progresses the phases rotate 
together with angular velocity to = (2iv/h)n s y%. The overlap 
between the orbits enforces the long-range phase coherence. 


a metal becoming superconducting expels the magnetic 
field from its interior. It is generally a very non-trivial 
question whether a many-body system will reach its low¬ 
est energy state predicted by thermodynamics. For ex¬ 
ample a ferromagnet cooled below its critical tempera¬ 
ture will not in general achieve an ordered state with 
macroscopic magnetization but rather break up into do¬ 
mains. A liquid when cooled will often end up as a solid 
glass or a polycrystal rather than a macroscopically or¬ 
dered single crystal. In these processes the forces at play 
are well studied and understood. In contrast it is re¬ 
markable that it is generally assumed that this question 
is not relevant to superconductors. How superconduc¬ 
tors achieve the state with the magnetic field expelled 
is not regarded to be an open question in the field [1]. 
The question has not even been posed, let alone been 
answered, in the vast literature dealing with the conven¬ 
tional theory of superconductivity. We suggest that the 
question has not been posed nor answered because the 
conventional theory lacks essential physical ingredients 
that are necessary to answer it. 

It has been argued that the explanation (or lack 
thereof) of the Meissner effect is in essence the same 
as that of flux quantization in superconducting rings 
[54, 55]. We argue that this is not quite so, even though 
the phenomena are certainly closely related. Flux quanti¬ 


zation involves changes in a fraction of the flux quantum 


with a = e 2 /(hc ) ~ 1/137 the fine structure constant. 
Eq. (64) is a microscopic quantity (~ 430 electron 
charges). How superconducting rings manage to adjust 
their current flow in order to respect flux quantization 
is certainly a fundamental question, as discussed exten¬ 
sively by A.V. Nikulov [54], for which we don’t have a dy¬ 
namical explanation. However, expulsion of a 200 Gauss 
magnetic field from a sample of cross section lew 2 is a 
much bigger question: it involves getting rid of some 10 9 
flux quanta! This is a macroscopic phenomenon for which 
it is reasonable to expect an explanation that is consis¬ 
tent with the macroscopic laws of physics. The fact that 
we don’t have a dynamical explanation of flux quanti¬ 
zation is not in our view a valid reason to argue that a 
dynamical explanation of the Meissner effect is not re¬ 
quired from the 60-year old theory that claims to explain 
conventional superconductivity [55]. 

It is very suggestive that the process of expulsion 
of magnetic field and transition to the superconduct¬ 
ing state of a macroscopic sample occurs through radial 
expansion of small superconducting regions (Fig. 1(b) 
or 1(c)), rather than uniformly as depicted in Fig. 1(a). 
We argue that this known experimental fact gives a 
vivid image of the underlying physics, and have provided 
here two complementarly closely related physical expla¬ 
nations: (i) expansion of a perfectly conducting fluid 
(Fig. 2(b)) leads to expulsion of magnetic field, as is 
well known in plasma physics (Alfven’s theorem [56, 57]), 
and (ii) expansion of electronic orbits leads to increased 
Larmor diamagnetism (Eq. (40)), as is well known in 
atomic and solid state physics [18]. In both processes the 
dynamical explanation for the magnetic field expulsion 
is the magnetic Lorentz force [58] acting on radially out¬ 
going charge. The driving force for a radial expansion 
is naturally a pressure, in our interpretation quantum 
pressure driven by reduction of quantum kinetic energy 
[38, 47, 59]. The postulated outflow of superconducting 
carriers and backflow of normal carriers (Fig. 4) strongly 
resembles processes known to exist in 4 He , where super¬ 
fluid thrusts from colder to warmer regions and normal 
fluid backflows from warmer to colder regions (fountain 
effect) [60]. In superfluid 4 He the driving force for the 
transition is known to be lowering of quantum kinetic en¬ 
ergy [59]. Finally, within our model the radial expansion 
of the phase boundary also gives a macroscopic image of 
the microscopic atomic orbital expansion that gives rise 
to superconductivity as described by the dynamic Hub¬ 
bard model [48]. 

In contrast, within conventional BCS theory there is 
no radial charge flow associated with the radial expan¬ 
sion of the phase boundary, hence the driving force for 
the azimuthal current cannot be the magnetic Lorentz 
force and remains unidentified. The driving force for the 
radial expansion of the phase boundary is termed ‘Meiss- 
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ner pressure’ by F. London [35] but is not given a physi¬ 
cal interpretation, and as a consequence the expansion of 
the phase boundary bears no relationship to the under¬ 
lying microscopic physics believed to be responsible for 
superconductivity within BCS [6], namely the Frohlich 
electron-phonon interaction and Cooper pairing. It could 
be said that conventional BCS theory describes the ‘ex¬ 
pansion’ of phase coherence as the superconducting nu¬ 
cleus expands. Then, BCS theorists have to explain how 
the expansion of phase coherence causes an azimuthal 
force to act on charge carriers. Perhaps this is a new 
force of nature that has not yet been identified [61] and 
will make BCS become part of the ‘reductionist frontier’ 
[62]. In our model instead, no new forces are needed 
and the development of phase coherence is explained by 
the expansion of the orbits to become strongly overlap¬ 
ping (Fig. 10), which in turn is directly linked to the 
azimuthal force acting on charge carriers (Fig. 9). 

We argue that since the conventional theory does not 
describe charge expulsion it cannot describe the Meiss¬ 
ner effect. More generally, a theory of superconductiv¬ 
ity that can explain the Meissner effect by expulsion of 
charge has to know the difference between positive and 
negative charge, just as superconductors do [49]. We ar¬ 
gue that this rules out all theories of superconductivity 
that are electron-hole symmetric, as most theories includ¬ 
ing the conventional theory are. Furthermore, since the 
expulsion of charge carries along an increase in poten¬ 
tial (Coulomb) energy, we argue that this rules out any 
theory of superconductivity where the condensation en¬ 
ergy is potential rather than kinetic which is the case of 
most theories including the conventional theory. Within 
our theory, electron-hole asymmetry and kinetic energy 
lowering are inextricably linked [48]. 

We have shown that the magnetization that the super¬ 
conductor develops to cancel the applied magnetic field 
originates in the orbital magnetic moments of electrons 
residing in mesoscopic orbits of radius 2 Al. The charge 
expulsion discussed in the previous paragraphs originates 
in the orbit enlargement from microscopic radius kp 1 2 3 4 5 to 
radius 2 Al in the transition to superconductivity [46]. 
In the absence of an applied magnetic field, electrons in 
these orbits give rise to a Spin Meissner effect [39], the 
existence of a macroscopic spin current within a Lon¬ 


don penetration depth of the surface of superconductors. 
The dynamical generation of the spin current and the 
expulsion of external magnetic fields occur through the 
same process, the expansion of the electron orbit from 
microscopic radius to the mesoscopic radius 2 Al in the 
presence of external magnetic field and internal electric 
field from the background ionic charge distribution. The 
fact that the speed of electrons in these orbits gives rise 
to angular momentum of value precisely h/2 [39] we re¬ 
gard as compelling evidence that the theory applies to 
nature. 

Finally we have pointed out that to understand both 
the Meissner and the Spin Meissner effects it is necessary 
to assume that different parts of the macroscopic super¬ 
fluid wavefunction 'F (r ) do not influence each other [43] , 
i.e. the magnetic field generated by electrons in overlap¬ 
ping orbits does not affect the magnetic field sensed by 
a given electron, and the ionic background electric field 
giving rise to the spin-Meissner current is not screened 
by the charge of electrons in overlapping orbits. This is 
consistent with the fact that the macroscopic wavefunc- 
tion of the superconductor T (r) is in many ways similar 
to the wavefunction 'F(r) of a single electron [26]. 

The process of negative charge expulsion and existence 
of mesoscopic orbits described here gives rise to a small 
charge inhomogeneity over the entire macroscopic sam¬ 
ple with the region within Al of the surface having a 
small excess negative charge. The resulting macroscopic 
equilibrium electrodynamics equations giving the spatial 
distribution and quantitative values of the charge den¬ 
sity, electric field and spin current in the ground state of 
superconductors are given in other publications [45, 46]. 
A valid microscopic theory of superconductivity will have 
to be consistent with these macroscopic and mesoscopic 
properties. 
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